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Abstract. In a companion paper, Maiolino et al. (2000) 
presented various observational evidences for "anomalous" 
dust properties in the circumnuclear region of AGNs and, 
in particular, the reduced Eb-v/Nh and Ay/Nn ratios, 
the absence of the silicate absorption feature in mid-IR 
spectra of Sy2s and the absence of the carbon dip in UV 
spectra of reddened Syls. In this paper we discuss various 
explanations for these facts. 

The observational constraints favor a scenario where 
coagulation, catalyzed by the high densities in the cir- 
cumnuclear region, yields to the formation of large grains. 
The resulting extinction curve is featureless, flatter than 
Galactic and the Eb-v/Nh and Ay/Nn ratios are signif- 
icantly reduced. These results should warn about an un- 
appropriate use of the standard Galactic extinction curve 
and Ay /Nh ratio when dealing with the extreme gas con- 
ditions typical of the circumnuclear clouds of AGNs. 

We also investigated alternative scenarios for the ob- 
served anomalous properties of dust in AGNs. Some of 
these scenarios might explain some of the observed prop- 
erties for a few objects, but they generally fail to account 
for all of the observational constraints obtained for the 
large sample of AGNs studied in these works. 

Key words: Galaxies: Seyfert - Galaxies: nuclei - Galax- 
ies: ISM - dust, extinction 



1. Introduction 

Some authors in the past have found evidences for anoma- 
lous properties of the dust in the circumnuclear region of 
AGNs. In particular, dust absorption and reddening ap- 
pear to be systematically lower than expected from the 
gaseous column density, for a Galactic gas-to-dust ratio 
and extinction curve (eg. Maccacaro et al. 1982, Reichert 
et al. 1985, Granato et al. 1997). In a companion paper 
(Maiolino et al. 2000, paper I) we have strengthened these 
previous results and provided additional evidences that 
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the dust in the circumnuclear region of AGNs is character- 
ized by anomalous properties. By comparing the redden- 
ing toward the Broad Line Region (BLR) and the gaseous 
N# inferred from the X-rays in a sample of 20 AGNs we 
found that, with the exception of three Low Luminosity 
AGNs, the Eb-v/Nh ratio is lower than the Galactic stan- 
dard value by a factor ranging from ^3 to ^100. By com- 
paring the optical and X-ray properties of various AGN 
classes (type 1.8-1.9 Sys, hard X-ray selected and radio 
selected QSOs, Broad Absorption Line QSOs and grism- 
selected QSOs) it was inferred that the Ay/Nn ratio must 
be one or two orders of magnitude lower than Galactic. 
Despite the large mid-IR continuum absorption inferred 
from the PAHs equivalent width, the absence of a signif- 
icant silicate absorption feature at 9.7/xm in the average 
ISO spectrum of Sy2s indicates that dust grains in the 
circumnuclear region of AGNs have different properties 
with respect to the Galactic diffuse interstellar medium. 
Finally, the lack of the 2 175 A carbon dip in the UV spec- 
tra of some type 1 AGNs, which appear affected by dust 
reddening, suggests a depletion of small grains in the ab- 
sorbing medium, as mentioned in paper I. These findings 
suggest that the dust grain distribution is biased in favor 
of large grains. The presence of large grains in the cir- 
cumnuclear region of AGNs was also suggested by Laor & 
Draine (1993) to explain some of the infrared properties 
of AGNs. 

In this paper we examine some possible explanations 
of the observational findings presented in paper I. In 
particular, we compare our models with the quantita- 
tive constraints given on the Eb_v/Njj ratio, the equiv- 
alent width of the silicate feature and the depth of the 
2175A carbon dip. Note that we often distinguish the 
Eb-v/Nh and Ay/Nn ratios since, when discussing ex- 
tinction curves different from Galactic, reddening and ab- 
sorption are no longer tied by the standard Galactic rela- 
tion Av/Eb-v = 3.1. In other words, models which could 
account for a decreasing of Eb-v/Nh do not necessarily 
lower the Ay/Nn ratio, and vice versa. 

In Section 2 we first analyze and, in practice, dismiss 
scenarios with a standard (Galactic) grain size distribu- 
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tion. The case for large grains is analyzed and discussed 
in Sect. 3 while in Sect. 4 we draw our conclusions. 

2. The case for standard grains 

2.1. Metallicity effects 

Various studies have found evidence for super-solar metal- 
licitics in the nuclear region of AGNs (Ferland et al. 
1996, Hamann & Ferland 1993, 1999, Fosbury et al. 1999, 
Reynolds et al. 1995, Dietrich & Wilhelm-Erkens 2000). 

The gaseous column densities used in paper I were in- 
ferred from the photoelectric cutoff observed in the X-rays, 
but actually such a cutoff gives a measure of the column 
of metals (Morrison & McCammon 1983), which is 
then converted into hydrogen column N# by assuming 
solar abundances. If the metallicity is higher than solar, 
then the N# inferred from the X-rays would be overesti- 
mated. On the other hand, a large fraction of refractory 
elements is locked into grains. As a consequence, an in- 
creased metallicity would also imply a higher dust con- 
tent. Both the dust content and the Nz to N# conversion 
factor depend linearly on the metallicity and, therefore, 
any change in metallicity should not affect the Eb-v/Nh 
ratio in the way we measure it. Finally, we should empha- 
size that there are no obvious astrophysical nucleosynthe- 
sis processes that favor the formation of non-refractory 
elements over refractory ones, implying that even "pecu- 
liar" metallicities (in terms of relative abundances between 
metals) should not affect the observed Eb-v/Nh ratio. 

Also, metallicity effects cannot explain the discrepancy 
between the lack of silicate feature and the strong mid-IR 
absorption in Sy2s. In the Small Magellanic Cloud the lack 
of the 2175A carbon dip is ascribed to a low metallicity, 
but this is unlikely to be the explanation for AGNs since, 
as discussed above, generally the latter are characterized 
by super-solar metallicities. 

Although metallicity effects cannot explain the anoma- 
lous properties of dust in AGNs, it can have important 
consequences on the expected emission line fluxes in the 
case of ionized gas containing dust. This is an issue which 
will be discussed further in sections 2.6 and 3. 

2.2. Reflected broad lines 

Another possibility is that the broad lines used to mea- 
sure Eb-v m paper I are not observed through the same 
medium that absorbs the X-rays (the obscuring torus), 
but are scattered by a medium (mirror) observed through 
a column density much lower than along the line of sight 
of the X-ray source. In this case the nuclei of these ob- 
jects are expected to show significant polarization. Tab. 1 
gives the optical and/or near-IR polarization for the ob- 
jects examined in paper I (therein Tab. 1) for which this 
information is available. For most of them the polarization 
is <2%, at least in the wavebands where the reddening is 
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0.6 


transm. 


2,9 



Table 1. Percentage of polarization in the optical (B,V) 
and in the near-IR (H) for the subsample of the sources 
listed in Tab. 1 of paper I for which this information 
is available. Polarization fraction in parenthesis arc cor- 
rected for dilution of the stellar continuum. Notes: a 
transm - objects for which there is substantial evidence 
that most of the polarization is due to dichroic transmis- 
sion; Pna~ gives the polarization percentage of the broad 
Ha line. References: 1) Barth et al. (1999), 2) Young et 
al. (1996), 3) Kay (1994), 4) Brindle et al. (1990), 5) Ber- 
riman et al. (1989), 6) Goodrich (1992), 7) Corbctt et al. 
(1998), 8) Wills et al. (1992), 9) Goodrich & Miller (1994). 



estimated (or as directly indicated by the polarization of 
Ha in two cases). For NGC5506, NGC2992 and Mkn231, 
which show higher polarization, there are evidences that 
most of the polarization is not due to reflection but to 
dichroic transmission. 

Also, the scattering efficiency generally does not ex- 
ceed a few percent. As a consequence, we would expect 
that the broad lines observed in our objects, once cor- 
rected for the observed reddening, are underluminous by 
a factor of ^100 with respect to the intrinsic luminosity 
indicated by the absorption-corrected X-ray luminosity, 
when compared to unobscured type 1 Seyfcrts and QSOs. 
Fig. 1 shows the broad Ha luminosity, reddening-corrected 
assuming a Galactic standard extinction curve, versus the 
intrinsic X-ray luminosity, again for the objects listed in 
Tab. 1 of paper I. The solid line and the shaded region 
show the \jB_a vs. hx relation (mean and ±ltr dispersion 
respectively) derived for a large sample of PG QSOs and 
nearby Syls not showing evidence for absorption. Clearly 
none of the objects is underluminous in broad Ha by the 
factor of ~ 100 expected in the case of reflection. The most 
"#a-deficient" object is NGC2992. This object does show 
significant polarization (Tab. 1), but as discussed above 
this is very likely a dichroic effect due to transmission 
through the host galaxy, which is seen edge-on; an inter- 
pretation which is also supported by the finding that the 
[OIII] line is polarized at the same level as the continuum 
and broad Ha line (Goodrich 1992). As we shall discuss 
in Sect. 3, the most likely interpretation for the Ha de- 
ficiency in this object is a peculiar extinction curve. It 
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Fig. 1. Ha luminosity for the objects reported in Tab. 1 of 
paper I, corrected for the reddening derived by assuming 
a Galactic extinction curve, as a function of the 2-10 keV 
luminosity (also corrected for absorption). The thin line 
and shaded region show the relation of these two quantities 
observed in Syl and QSOs (mean and ±ler dispersion). 



should also be noted that Mkn231 is actually significantly 
over luminous in Ha, by a factor of ^30. This object is 
actually known to be underluminous in the X-rays. Possi- 
bly, the absorbing column density might be much higher 
than that estimated by Iwasawa (1999) and Turner (1999) 
in their faint hard X-ray spectrum. In the latter case the 
intrinsic X-ray luminosity could be much higher. 

We should however mention that the sample presented 
in paper I (Tab. 1) is probably biased in favor of objects 
which are particularly bright in the broad lines emission, 
just as a consequence of the selection criterion (i.e. those 
absorbed objects which show at least some broad emis- 
sion lines). Nonetheless, it is highly unlikely that such a 
selection effect accounts for a factor of ~100 in the Ha 
flux. 

It is worth noting that a diagram analogous to Fig. 1 
was presented by Veilleux et al. (1999, therein Fig. 3) when 
discussing the detection of broad Paa in Ultraluminous In- 
frared Galaxies (ULIRGs); there the broad Ha is replaced 
by the broad Paa and the (intrisic) X-ray luminosity is 
replaced by the far-IR luminosity. From that diagram it 
can be similarly inferred that the broad component of the 
Paa detected in ULIRGs cannot be ascribed to reflection. 
For what concern Seyfert galaxies, Veilleux et al. (1997) 
reach a similar conclusion stating that the broad compo- 
nent of the hydrogen near-IR lines detected in some Sy2s 
cannot be ascribed to reflection. 



Finally, also in the reflection scenario the lack of sili- 
cate absorption in Sy2s and of the lack of the carbon dip 
in reddened Syls is not explained. 

2.3. Reduction of the extinction scattering component 

Unlike the case of Galactic stars where both absorption 
and scattering contribute to the extinction, in obscured 
AGNs we may have effects of scattering back into our 
beam. This would imply a reduction of the scattering com- 
ponent in the effective extinction curve. However, the net 
effect must be small compared to what observed for the ob- 
jects of our sample. Even in the most extreme case that the 
absorbing/scattering medium covers the whole solid angle 
seeing by the nuclear source (that is very unlikely) and 
that the scattered radiation enters our beam without un- 
dergoing any dust absorption (that is also very unlikely), 
then the effect would be to completely remove the scatter- 
ing component out of the total exinction curve. Since at 
the wavelengths of interest T scat ~ 2 r a b s ~ 2/3 r ext , the 
effect of this extreme case would be to reduce the extinc- 
tion and reddening by a factor of 3 at most. While in our 
sample we observe several cases for which the reddening 
deficit is much higher. 

The analysis above makes the implicit assumption that 
the line of sight is "typical", but it is possible that the 
absorbing column along the line of sight is lower than the 
scattering column along the other directions. In the latter 
case, and under the extreme assumption discussed above, 
the effective extinction could be reduced further. Yet, this 
effect of scattering back into the beam would not explain 
the lack of absorption features in the UV and in the mid- 
IR. 

24. Extended BLR 

Several reverberation mapping studies have shown that 
the BLR has a size of ~0.01 pc. More specifically R ss 

1 /2 

0.02 L 4 g pc (Clavel et al. 1991, Peterson 1993, Salamanca 
et al. 1994). On the other hand, various pieces of evidence 
indicate that the obscuring torus has a size of ~1 pc or 
larger (Gallimore et al. 1997, Greenhill & Gwinn 1997, 
Ford et al. 1997, Granato et al. 1997, Maiolino 2000, and 
references therein) . Since the scale of the latter is so much 
larger than both the X-ray source and the BLR, the ab- 
sorbing medium along the line of sight should be the same 
for both. However, for many of the objects in the sample 
presented in paper I (Tab. 1) the broad lines have a width 
of about 2000 km/s, which is somewhat narrower than 
usually observed in type 1 Seyferts and QSOs. This might 
imply that in these objects the BLR, or the fraction of it 
that we observe, is more extended than assumed and with 
dimensions comparable to the obscuring torus. As a conse- 
quence, the BLR might be observed, on average, through a 
lower absorbing column with respect to the nuclear X-ray 
source. 
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However, since the broad lines in many of these objects 
are weak, very broad wings might have been undetected 
because either lost in the noise or confused with the un- 
derlying continuum features, thus yielding to an under- 
estimate of the real width of these lines. Also, in some of 
these objects the broad lines were observed to vary on time 
scales shorter than a few months, suggesting a BLR scale 
smaller than 0.1 pc (Eracleous & Halpern 1993, Winkler 
et al. 1992, Rosenblatt et al. 1994, Giannuzzo & Stirpe 
1996, Peterson et al. 1984, 1982, Miller et al. 1985, 1979). 

Summarizing, we cannot exclude that for some of the 
objects in our sample the BLR is extended enough that the 
broad lines and the (much smaller) hard X-ray source are 
obscured by columns of gas that are significantly different. 
However, this is unlikely to be the case for most of the 
objects. 

Finally, the lack of the silicate feature in the mid-IR 
spectra of Sy2s and the lack of the carbon dip in the UV 
spectra of reddened Syls remain unexplained by the dif- 
ferential absorption scenario discussed in this section. 

2.5. Absorption by BLR clouds 

If one or more BLR clouds are located along the line of 
sight toward the X-ray nuclear source then they would 
contribute to the observed gaseous absorbing column den- 
sity (BLR clouds are expected to have Nh ~ 10 23 cm~ 2 ), 
but obviously this absorbing column would not affect 
the observed broad line flux. However, should the low 
Eb-v/Nh be a feature common to most AGNs and should 
the cause of this effect be absorption by BLR clouds, 
then this would require a large covering factor of the BLR 
clouds, at least as high as the covering factor of the ob- 
scuring medium [], i.e. ~80% . On the other hand the 
covering factor of the BLR clouds is estimated to be low 
(~10%) based on the absence of any Ly-edge cutoff in the 
UV spectrum of QSOs. Therefore, absorption due BLR 
clouds might explain some of the cases of low Eb-v/Nh, 
but probably not the majority of them. 

Another reason to discard this explanation is related 
to the lack of strong and short term variability of the 
absorbing N# in the X-rays. If most of the photoelectric 
cutoff observed in the hard X-rays is due to a single or 
a few BLR clouds along the line of sight, then a strong 
variation of the Nh would be expected on a time scale of 
a few weeks or less, given the expected size of the BLR 
clouds (< 0.01 pc) and their velocity. However, about 1/3 
of the objects in our sample were observed two or more 
times in the hard X-rays (eg. Gilli et al. 2000, Malizia et 
al. 1997) and the Nh inferred from the cutoff does not 
show the strong variations expected from the motion of a 
BLR cloud passing on the line of sight. There are a few 

1 A covering factor of 80% for the obscuring torus is inferred 
both from the Sy2/Syl ratio and from the opening angle of the 
light cones observed in Sy2s (Maiolino & Rieke 1995). 



objects (most of them not in our sample) whose hard X- 
ray monitoring has indeed shown variations of Nh, but 
such variations are on time scales of years, i.e. consistent 
with a distance of the absorber from the nucleus of a few 
pc (Maiolino 2000). The only case known to show rapid 
N H variability is NGC4151. 

Again, the issues related to the silicate feature and to 
the carbon dip remain unexplained also in this case. 

2.6. Dust sublimation 

The inner part of the obscuring torus is directly exposed 
to the UV radiation emitted by the nuclear source. If the 
UV flux at this location is large enough, the dust equilib- 
rium temperature can reach the sublimation value (Laor 
& Draine 1993, Phinney 1989). If dust is sublimated for a 
significant fraction of the absorbing gaseous column this 
would explain the reduced Eb-v with respect to that ex- 
pected from the observed N#. However, this model has 
problems when compared with the low covering factor in- 
ferred for the ionized gas in the circumnuclear region of 
AGNs, as discussed in the following. 

The capability of the dust in competing with the gas 
in the absorption of UV ionizing photons increases signif- 
icantly with the ionization parameter U. Indeed, a highly 
ionized gas is much more transparent than dust to the UV 
ionizing photons. In particular, for a Galactic gas-to-dust 
ratio and dust composition, dust grains are expected to 
dominate the absorption of photons beyond the Lyman 
edge when U is higher than about 10~ 2 . At the inner face 
of the torus the ionization parameter is expected to be 
about 0.1 or larger (Netzer & Laor 1993, Netzer 2000) 
and therefore most of the UV ionizing photons should 
be absorbed by the dust which reprocesses this radiation 
into near and mid-IR emission, which is actually observed 
in the nuclear region of AGNs (eg. Maiolino et al. 1995, 
Clavel et al. 2000, Maiolino et al. 1998, Thatte et al. 1997, 
Oliva et al. 1999). As a consequence, although the obscur- 
ing torus probably subtends as much as 80% of the solid 
angle seen by the nuclear UV source, very little line emis- 
sion is produced. The emission lines that would be pro- 
duced by the torus in absence of dust should have a width 
significantly smaller than the broad lines and close to the 
width of the narrow lines. The effect of strong absorption 
of UV ionizing photons from dust at the inner face of the 
torus has been proposed to explain the marked separation 
between BLR and NLR and to explain the low covering 
factor of the NLR. A more detailed discussion of this is- 
sue is given in Netzer & Laor (1993), in Laor & Draine 
(1993) and Pier & Voit (1995). Here we do not discuss 
further the implications of dust absorption for the pho- 
toionization models of the circumnuclear gas, but we use 
these considerations as a constraint for the properties of 
the dust at the inner face of the torus. 

If most of the gas in the torus is dust-free, especially at 
the inner face, because of dust sublimation, then photons 



R. Maiolino, A. Marconi, & E. Oliva: Dust in active nuclei 



5 



at wavelengths shortward of the Ly edge should ionize the 
gas in the torus and should produce powerful hydrogen 
lines nearly as narrow as the narrow lines. The covering 
factor of the torus is estimated to be ~80%. Until a few 
years ago the covering factor of the NLR was estimated to 
be ~2% (eg. Netzer & Laor 1993); however recent obser- 
vations of the UV spectrum of QSOs indicate a shortage of 
ionizing photons with respect to former estimates by a fac- 
tor of about 4 (Zheng et al. 1997, Laor et al. 1997) there- 
fore implying a larger covering factor of the NLR (^8%), 
but which is still much lower than the covering factor es- 
timated for the obscuring torus. As a consequence, the 
flux of the emission lines produced by the dust-sublimated 
torus should be about 10 times larger than the observed 
narrow lines. 

A possible way to suppress the line emission in the 
dust sublimated medium is to make the ionized gas very 
hot. However, in the latter case the absorbing medium 
observed in the hard X-rays should appear "warm" (ac- 
cording to the X-ray standards). On the other hand, in 
most of the objects of the sample presented in paper I the 
absorbing gas is "cold" beyond doubt. 

Summarizing, it is difficult to ascribe the low 
Eb-v/Nh ratio to dust sublimation at the inner face of the 
torus without running into a serious emission line budget 
problem. 

For what concerns the absorption features, dust subli- 
mation might explain the lack of the carbon dip at 2 175 A 
in reddened Syls, since the very small grains which are re- 
sponsible for this feature undergo sublimation more eas- 
ily than larger grains. However, it is difficult to ascribe 
the lack of the silicate feature in the spectra of Sy2s to 
sublimation, indeed the absence of this absorption feature 
requires a dust distribution dominated by grains larger 
than a few /im (i.e. significantly larger than the maximum 
dimension of 0.25/mi inferred for the Galactic medium), 
rather than the destruction of small grains (Sect. 3. 3). 

2. 7. Low dust-to-gas ratio 

One of the simplest ways to explain the observed low 
Eb-v/Nh is to assume that the dust-to-gas mass ratio in 
the obscuring medium is lower than the standard Galac- 
tic value. However, a reduced abundance of dust in the 
dense medium of the circumnuclear clouds is not obvious 
to justify. The vicinity of the AGN could help to create the 
conditions that act to destroy dust grains. In particular, 
sputtering by X-rays and/or shocks might be a possible 
mechanisms of dust destruction in this environment. 

Obviously, the Eb-v/Nh ratio decreases linearly with 
the dust-to-gas ratio, if the dust properties remain un- 
changed. This is shown in the top panel of Fig. 2. The 
trend of Eb-v/Nh as a function of the dust-to-gas ratio 
relative to the Galactic value is given by the thick solid 
line. The horizontal thin line and shaded region show the 
mean and ±lcr dispersion of Eb-v/Nh for the objects in 



the sample presented in paper If]. A factor of about 10 
lower dust-to-gas ratio is required to match the mean of 
the observed Eb-v/Nh values. The thin dashed vertical 
lines indicate the range of gas-to-dust ratios that matches 
the observed Eb-v/Nh for ^68% of the objects (±lc), 
which implies a dust-to-gas ratio between ~ 3 times and 
50 times lower than Galactic. A dust underabundance of 
factor of ^100 is required to explain the lowest Eb-v/Nh 
found in paper I (see Fig. 1 in that paper). The central 
panel shows the trend for Ay/Nn, which also decreases 
linearly since the absorption curve was assumed constant. 

As the dust content decreases, the absorption of UV 
ionizing photons by dust grains becomes gradually less 
important and, as discussed in the previous section, the 
emission of nebular lines gradually increases. In order to 
check whether the required dust depletion runs into the 
same emission line budget problem as in the case of the 
sublimated dust (Sect. 2. 5), we have created a simple semi- 
-analytical code to calculate the ionization structure and 
hydrogen lines emissivity of a gaseous cloud containing a 
variable amount of dust. We assumed a ionization param- 
eter U=0.1 at the inner face of the cloud, as this is about 
the value expected at the inner face of the torus (Netzer & 
Laor 1993). The covering factor was assumed to be 80%. 

The results of the model were then compared with 
the expected emission line flux in the case of a dust- 
free medium with a covering factor equal to unity, i.e. 
the maximum hydrogen line flux (H/3(max)) which can be 
obtained with the available photoionizing luminosity. The 
ratio between the hydrogen line flux in the dusty and in 
the dust-free cases gives the importance of the line emis- 
sion from the obscuring torus relative to other gaseous 
clouds which are either free from dust, such as the BLR 
clouds, or whose ionization parameter is low enough to 
make unimportant the effect of dust absorption of UV ion- 
izing photons, such as the NLR clouds. The thick black 
solid line in the bottom panel of Fig. 2 shows the ratio 
between the H/3 flux in the case of a dusty medium (80% 
covering factor) and the case of a dust-free medium (100% 
covering factor), H/3(torus)/H/3(max), as a function of the 
dust-to-gas ratio in the obscuring torus. 

We also performed the same calculation with Cloudy, 
the photoionization code provided by Ferland and collab- 
orators (Ferland 1999), which accounts for several effects 
which are not included in our simple code. The output of 
Cloudy is overplotted with a thin solid line in the bottom 
panel of Fig. 2. Cloudy's result is not much different from 
ours and, partly, it is also to ascribe to the extinction curve 
used by Cloudy which is somewhat different from the 
Galactic standard curve (eg. Ry = Av/Eb-v = 4, to be 
compared with the standard Galactic value of Ry = 3.1). 
This comparison justifies us in using our code also for 

2 In this analysis we excluded the three Low Luminosity 
AGNs and the three objects whose reddening was estimated 
by means of the continuum fit (paper I). 
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other models in the next section (Sect. 3). The advan- 
tage of our code is that is it possible to use an arbitrary 
dust composition with a (quasi-) continuous distribution 
of grain sizes. 

The H/3(torus)/H/3(max) values obtained with the 
model discussed above have to be compared with the 
flux of the narrow lines. More specifically, the flux of 
the lines emitted by the torus should be lower than the 
flux of the narrow lines not to incur in the emission 
line budget problems discussed in Sect. 2. 5; this implies 
that H/3(torus)/H/3(max) must be lower than the cov- 
ering factor of the NLR clouds, i.e. <0.08 (Sect. 2. 5). 
This limit is shown with an horizontal line in the bot- 
tom panel of Fig. 2. The dust-to-gas ratio required 
to match the mean observed Eb-v/Nh would give 
H/3(torus)/H/3(max)«0.15 that is inconsistent with the 
NLR covering factor. Although, higher values of Eb_ v /Nh 
might be marginally consistent with the NLR emission, 
the values of Eb-y/Nh lower than the mean would imply 
an even higher H/3(torus)/H/3(max), as shown in Fig. 2. 

In Sect. 2.1 we mentioned that the nuclear region of 
AGNs is generally characterized by super-solar metallici- 
ties. As we are going to discuss, a high metallicity helps 
to relax the constraints on the emission line flux. Metal- 
licity variations change the dust-to-gas ratio, but the vari- 
ation of the dust-to-gas ratio due to metallicity leaves 
unaffected the measured Eb-v/Nh and Ay/Nn- As dis- 
cussed in Sect. 2.1, this is because the Nh is measured in 
the hard X-rays through the column of metals. On the 
other hand, the increased dust-to-gas ratio, regardless of 
whether due to metallicity or to other effects, does af- 
fect the dust-to-gas opacity for the UV ionizing photons 
and, therefore, changes the H/3(torus)/H/3(max) accord- 
ingly. As a consequence, even if metallicity does not af- 
fect the Eb-v/Nh ratio, it can "decouple" the latter ratio 
from the emission line budget constraints. In particular, 
the dotted lines in Fig. 2 show the effect of a metallic- 
ity three times higher than solar. On the top and central 
panel of Fig. 2 the (metallicity-) increased dust-to-gas ra- 
tio leaves unaffected the quantities on the Y-axis, thus 
yielding a shift on the X-axis. Instead, in the bottom panel 
the H/3(torus)/H/3(max) ratio does change accordingly to 
the reduced dust-to-gas ratio, i.e. this ratio keep following 
the trend given by the solid line. As a result, for the same 
Eb-v/Nh ratio the implied H/3(torus)/H/3(max) is signif- 
icantly lower. In particular, at the dust-to-gas ratio which 
gives an Eb-v/Nh accounting for the mean observed 
value, H/3(torus)/H/3(max) as 0.07, which is marginally 
consistent with the covering factor of the NLR. Yet, for 
lower values of Eb-v/Nh the line emission from the torus 
is still much larger than the NLR. 

Summarizing, a low gas-to-dust (mass) density ratio 
might probably explain the reduced Eb-v/Nh for sev- 
eral of the objects without incurring in emission line bud- 
get problems, especially if the metallicity is super-solar 
as claimed by several authors. However, low values of 



Eb-v/Nh (lower than the mean) are still difficult to ex- 
plain with this model. 

Certainly, the reduced dust-to-gas ratio alone cannot 
account for the lack of the silicate feature in the mid- 
IR spectra of Sy2, nor for the lack of the carbon dip in 
reddened Syls. 

3. The case for large dust grains 

3.1. Plausibility of the large grains scenario 

So far we have adopted the standard extinction curve ob- 
served in the diffuse Galactic interstellar medium. How- 
ever, the physical conditions of the gas in the nuclear re- 
gion of Seyferts and QSOs are so extreme that there is 
no reason to assume that dust has the same properties 
as in the diffuse interstellar medium of our Galaxy. The 
gas densities of the obscuring medium in AGNs probably 
exceeds 10 5 -10 6 cm~ 3 ,[] i.e. a density similar or higher 
than observed in the densest dark clouds of our Galaxy. 
Along the line of sight of dense Galactic molecular clouds 
the extinction curve is generally found to be flatter than 
in the diffuse ISM (Vrba & Rydgren 1984, Carrasco et 
al. 1973, Whittet et al. 1976, Cardelli et al. 1989). Vari- 
ous observational evidences indicate that most likely such 
flat extinction curves are due to a distribution of grain 
sizes on average larger than in the diffuse ISM (Kim et 
al. 1994, Weingartner & Draine 2000). Larger grains may 
form either by accretion from the gas phase or by grain 
coagulation. The former case should increase the Ay/Nn 
ratio, while the latter case should generally yield a reduced 
Ay/Nn- Although measurements of Ay/Nn are very un- 
certain, some works suggest that Ay/Nn is lower in dense 
clouds characterized by flatter extinction curves, there- 
fore supporting the coagulation scenario (Jura 1980, Kim 
& Martin 1996). Also, theoretical considerations indicate 
that coagulation is much more effective for the growth 
of dust grains than accretion from the gas phase (Draine 
2000). 

Since the grain coagulation rate increases with density 
(oc n 1 / 2 , Draine 1985), the effect is expected to be even 
more dramatic in the circumnuclear clouds of AGNs. In 
particular, for a gas density in excess of 10 6 cm -2 the es- 
timated time scale for the depletion of small grains due to 
coagulation is shorter than about 10 5 yr, i.e. shorter than 
the typical dynamical time scale in the circumnuclear re- 
gion of AGNs. This not only implies that dust grains have 
time to coagulate before being replaced by gas coming 

3 Such a high density is inferred by the requirement of con- 
fining columns of gas higher than 10 24 — 10 25 cm -2 in a region 
smaller than a few pc, to avoid a gas mass higher than the 
dynamical mass (Risaliti et al. 1999). The presence of clouds 
with densities of 10 5 — 10 6 cm~ 2 in the circumnuclear region of 
AGNs is also directly inferred from HCN and other millimetric 
observations of AGNs with arcsecond resolution (eg. Tacconi 
et al. 1994). 
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from the outer regions, but actually some mechanism that 
prevents grains to grow without limits must be invoked 
(eg. shocks). 

Summarizing, a dust distribution biased in favor of 
large grains, probably as a consequence of coagulation, is 
naturally expected in the extreme physical conditions of 
the gas in the circumnuclear region of AGNs. 

The presence of large grains in the obscuring medium 
of AGNs was also proposed by Laor & Draine (1993) in re- 
lation to the IR properties of active nuclei. In the following 
we discuss the effects of large grains, specifically formed 
by coagulation, on the observational quantities discussed 
in paper I. 

3.2. The effect of large grains on Eb-v/Nh 

To obtain a more quantitative description of the effects 
of large grains in the circumnuclear region of AGNs we 
modelled the size grain distribution with a power-law: 

dn(a) oc a~ q da, a min < a < a max (1) 

where a is the grain diameter (assumed spherical) and 
dn(a) the density of grains with size between a and a+da. 
Absorption and scattering cross sections were calculated 
by using the refractive indices in Draine & Lee (1984) 
and in Laor & Draine (1994), and by numerically solving 
the Mie formulas. The extinction curve in the Galactic 
diffuse ISM was successfully described by Mathis et al. 
(1977) with a mixture of graphite (47%) and silicate (53%) 
grains and with a m i n = 0.005/im, a max = 0.25/im and 
q=3.5 0. The effect of large grains was investigated by 
us in various ways: either by increasing a max , or a m i n , or 
by decreasing q, or by a combination of these possibilities. 
Coagulation was assumed by imposing the total dust mass 
to be constant (relative to the gas mass). The increased 
average grain size always gives a flatter extinction curve. 
As long as a max is larger than 0.5^m, coagulation also 
yields to Eb-v/Nh and Ay/Nn lower than the Galactic 
standard value. We cannot discuss here all of the possible 
cases which give flatter extinction curves. In the following 
we discuss three representative families of models, all of 
which assume the = 0.005/im, but different 

a max and q. 

3.2.1. q=3.5, increasing a max 

Fig. 3 is analogous to Fig. 2, but in this case the dust- 
to-gas (mass) ratio is kept constant, while a max increases 
along the abscissa and the other two dust parameters are 
kept constant at q=3.5 and a m j n = 0.005. 

4 Note that this model does not include the very small grains 
(PAHs) which are thought to be responsible for many of the 
emission feature in the mid-IR spectrum and probably con- 
tribute significantly also to the absorption dip at 2175A (Wein- 
gartner & Draine 2000). 



In addition to the panels shown in Fig. 2, we also 
include in Fig. 3 a panel giving the variation of Rv = 
Av/Eb-v (the larger is Ry the flatter is the extinction 
curve). In this case Ry increases up to ~7 for very large 
a-max ■ Such large values of Ry imply that the absorption 
in the optical is higher than that estimated from the red- 
dening assuming a Galactic extinction curve. 

The Ha luminosities for the objects shown in Fig. 1 
were corrected assuming Ry=3.1. If Ry is higher than 
Galactic, then the Ha luminosities have been underesti- 
mated by the following quantity: 

ALog(L Hce ) - ^?=p (kR v - 3.1) 0.79 (2) 
2.5 

where 1 < k < 1.26 depending on the details of the ex- 
tinction curve. For an average Eb-v ~ 0.5 the correction 
is significant, i.e. w 0.7 in log L Ha . Some of the objects 
in Fig. 1 do appear to be that or even more underlumi- 
nous in Ha with respect to the Syl+QSO relation, thus 
supporting the large grains scenario. However, other ob- 
jects are not underluminous in Ha. For some of them 
the required L(Ha) correction is small and therefore still 
consistent with the Syl+QSO distribution, but for others 
the corrected Ha luminosity becomes uncomfortably high 
(~ 2a above the Syl+QSO distribution). However, as dis- 
cussed in Sect. 2. 2 our sample is probably biased in favor 
of Ha bright sources and, therefore, the excess of Ha lu- 
minosity discussed above does not necessary rule out the 
large grains hypothesis for these objects. 

The second panel of Fig. 3 shows the expected trend for 
Eb-v/Nh (thick solid line). This ratio has a steep decline 
as a function of a max . It should be noted that the locus 
of the values observed in AGNs (shaded area) is no longer 
constant as in Fig. 2, this is because we also took into 
account that the reddening estimate Eb-v based on the 
broad line ratios (paper I) changes with the extinction 
curve. The mean Eb-v/Nh of the AGNs in our sample 
can be explained with a max ~ 30/xm. A much larger value 
of a max (> 100/im) is required to match the properties of 
objects with lower E B _v/N H . 

The ratio Ay/Nn also decreases rapidly with a max , 
indicating that large grains can, more generally, explain 
the low Ay /Nh inferred for many AGNs and discussed in 
paper I. 

The emission line budget problem does not improve 
much with respect to the case of low dust-to-gas ra- 
tio discussed in Sect. 2. 6: for solar metallicities, the 
H/3(torus)/H/3(max) ratio ranges from ^0.04 up to > 0.25 
for the values of a max required to match the observed 
Eb-v/Nh. 

In the case of super-solar metallicities the dust-to-gas 
ratio increases, but both the extinction curve, and the 
Eb-v/Nh and Ay/Nn ratios remain unchanged, as dis- 
cussed in Sect. 2. 6. However, the H/3(torus)/H/3(max) ratio 
decreases as a consequence of the increased dust absorp- 
tion of UV ionizing photons (see discussion in Sect. 2.6). 
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The case of an increased dust-to-gas ratio due to a metal- 
licity three times higher than solar is shown with a dot- 
ted line in the bottom panel of Fig. 3. In this case the 
constraints given by the narrow line flux are met better 
and, in particular, for a max lower than about 100 fim, the 
expected H/3(torus)/H/3(max) is lower than the NLR con- 
straints. Yet, very low values of Eb_v/Nh (observed for 
some of the objects analyzed in paper I) remain difficult 
to explain with this model. 

3.2.2. a max = 1/zm, decreasing q 

In Fig. 4 we show another "family" of grain distribution 
models: here a max is kept constant and equal to 1/im, while 
q decreases along the abscissa. The extinction curve flat- 
tens very rapidly with q, as also indicated by the rapid 
growth of Ry. As shown, in the second panel, such flat 
extinction curve would explain the objects with nearly no 
reddening Eb-v, although substantial gaseous column Njj 
is observed along the line (paper I). However, several of 
the objects considered in paper I do show some redden- 
ing, therefore indicating that in these objects the extinc- 
tion curve cannot be completely flat. A possibility, is that 
in the outer parts of the obscuring torus, where the gas 
density is probably much lower, the properties of dust are 
more similar to those of the diffuse Galactic ISM. In this 
case the extinction, curve should be nearly identical to 
Galactic, but shifted to higher values of absorption as a 
consequence of the higher density part of the torus which 
contributes with a flat extinction curve. The presence of 
"normal" dust in the outer parts obviously reduces Ry, 
but also increases Eb-v/Nh- The case in which 15% of 
the gaseous absorbing column Nh contains "normal" dust 
is shown with a thick dashed line in Fig. 4. In this case 
Rv is not unreasonably high, and the Ay inferred from the 
observed reddening is still consistent with the presence of 
weak broad lines in the optical or IR spectra of the AGNs 
examined in paper I. Yet, the large value of Ry makes the 
Ha luminosity correction even higher than in Sect.3.2.1, 
which complies with the Ha-weak objects, but might be 
a problem for those objects which are above the distri- 
bution of Syl+QSO in Fig. 1 (see analogous discussion 
in Sect.3.2.1). Another problem with the inclusion of an 
even small (15%) fraction of dust with Galactic extinction 
curve is that the ratio Eb-v/Nh has a much slower de- 
cline and, as shown in Fig. 4, cases of very low Eb-v/Nh 
arc difficult to explain. 

Another important effect of this "family" of models is 
that, at variance with the Eb-v/Nh ratio, the Ay/Nn 
ratio is not reduced by a large factor, as shown in Fig. 4. 
Therefore, this class of models does not provide the most 
likely explanation for objects characterized by very low 
Ay/Nn values (paper I). 

Yet, one of the most appealing properties of this fam- 
ily of models is that the H/3(torus)/H/3(max) ratio is low 
(~8% at most) and may be much lower than the NLR 



constrain, especially if the dust-to-gas ratio is increased 
as a consequence of super-solar metallicity (bottom panel 
of Fig. 4). The value of H/3(torus)/H/3(max) is obviously 
unaffected by a fraction of the column density containing 
"normal" Galactic dust, if the latter is located in the outer 
parts of the torus. 

3.2.3. a max = lO^rn, q decreasing 

This family of models is shown in Fig. 5. These are similar 
to those discussed in Laor & Draine (1993). Here we do 
not discuss in detail the properties of these models, but 
mostly focus on the main differences and similarities with 
the former case. 

In these models the effect of decreasing q not only 
rapidly flattens the extinction curve, but also decreases 
significantly the absorption Ay, and therefore can easily 
explain little reddened objects with low Ay/Nn (paper I). 
If the metallicity is three times solar, or higher, the range 
of q required to match the Eb-v/Nh distribution of most 
objects (q>3.1) also complies with the NLR constraint on 
the torus emission line flux (bottom panel of Fig. 5). 

As a consequence of the much reduced absorption Ay, 
the inclusion of an outer layer of Galactic dust reduces 
much more Ry and affect much less Eb-v/Nh- This is 
shown by the thick dashed line in Fig. 5, which gives the 
case with only 3% of the (outer) column of gas containing 
Galactic "standard" dust. Therefore, this family of mod- 
els relaxes the problems on the large Ry discussed in the 
family of models discussed in Sect. 3. 2. 2: it gives very low 
Eb-v/Nh and also Ay/Nn values without incurring in 
large values of Ry . 

3.2.4. Summary of the effects on Eb-v/Nh 

We cannot describe here all of the dust models that 
matches the specific constraints for each single object. 
However, for most of the AGNs discussed in paper I one 
or more dust distribution models match the specific con- 
straints on Eb-v/Nh, Eb-v, broad lines luminosity and 
detection. 

3.3. The effect of large grains on the silicate feature 

As mentioned in the introduction and discussed more in 
detail in paper I, Sy2s appear characterized by significant 
mid-IR absorption which, nonetheless, is not associated 
to a deep silicate absorption feature at 9.7/um which is ex- 
pected by a Galactic extinction curve. Since this silicate 
feature is produced by grains smaller than a few /im, a 
dust distribution dominated by large grains can also ex- 
plain this "anomaly" of the mid-IR spectra of Sy2s, as 
discussed in the following. 

We take the optical depth at 7.7/im as representative 
of the featureless absorption in the vicinity of the silicate 
feature. We specifically choose this wavelength because is 
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the location of the PAH whose equivalent width is used 
by Clavel et al. (2000) to determine the mid-IR absorp- 
tion. The difference of the optical depths 79.7 — 77.7 is an 
alternative way to measure the equivalent width of the 
silicate feature, as long as it is not saturated. The ratio 
( T 9.7 — T i.7)/ T i.7 is a measure of the depth of the silicate 
feature relative to the featureless mid-IR absorption. In 
Fig. 6 the upper horizontal line indicates the standard 
Galactic value of this ratio, while the lower horizontal line 
indicates the (mean) upper limit of this ratio for Sy2s, 
as inferred from the data in Clavel et al. (2000), and the 
shaded area gives the spread of the upper limits (as in- 
ferred from the spread of the EW of the PAH at 7.7/xm). 
Fig. 6 illustrates in different terms what already discussed 
in paper I: the upper limit on the depth of the silicate 
feature in the average spectrum of Sy2s is about one or- 
der of magnitude lower than expected from the Galactic 
extinction curve. The curves plotted in the diagram show 
the expected (T9.7 — T7. 7)^7.7 ratio according to the co- 
agulation models discussed in Sect. 3. 2. In particular, the 
dashed line refers to the model presented in Sect. 3. 2.1 (q 
constant and equal to 3.5) and is related to the upper 
abscissa scale, while the solid curves refer to the models 
presented in Sects. 3. 2. 2 and 3.2.3 (a max constant) and are 
related to the lower abscissa scale. All of these models pre- 
dict a depth of the silicate feature significantly lower than 
Galactic. The model with a max = lOyum can easily fit the 
constrain given by the mid-IR observations, with q just 
slightly lower than the canonical value of 3.5. 

As mentioned in paper I (sect. 4.1), the inner (hotter) 
mid-IR emitting region probably is not featureless but is 
also characterized by the 9.7/im feature in emission. Since 
the upper limit on T9.7 — T7.7 was determined assuming an 
intrinsic featureless continuum, the possible presence of 
the 9.7/im feature in emission should relax such an upper 
limit. However, as mentioned in paper I, the equivalent 
width of the 9.7/xm emission feature is so small that the 
upper limit on (79.7 — j) / tj would increase by a quan- 
tity of 0.15 at most, which certainly cannot account for the 
the discrepancy with the value expected for the Galactic 
case (Fig. 6). 

Finally, we wish to emphasize that the presence of 
strong PAH features in the ISO spectra of Scyfcrt galax- 
ies is not in contractdiction with our large grain model, 
since such PAH features most likely come from the host 
galaxy (where the dust properties are probably similar to 
the Galactic ISM) and not from the nuclear region. 

3.4- The effect of large grains on the carbon dip 

In Sects. 3. 2 and 3.3 we only considered the effects of coag- 
ulation by increasing a max and by decreasing q, since these 
are the parameters which primarily affect Ry, Eb-v/Nh, 
Av/Nr and the silicate feature. The latter quantities are 
instead little sensitive to a m i n . However, this is not the 
case for the carbon dip which, in the model of Mathis 



et al. (1977) is ascribed to graphite grains with sizes in 
the range lOOA- 200A. In particular, the carbon dip is 
washed away from the extinction curve if coagulation de- 
pletes small grains by increasing a m i n above 200A. This is 
shown in Fig. 7, where we plot the variation in the UV 
shape extinction curve by modifying the dust distribution 
to have a m i n = 0.03/im and leaving a max and q to the 
standard values of 0.25 and 3.5, respectively. 

Note that recent studies have shown that most of the 
profile of the 2 175 A dip is probably produced by even 
smaller grains (PAHs, Weingartner & Draine 2000) that 
are not included in our model. Obviously, in this case the 
depletion of small grains onto larger grains through co- 
agulation would be even more effective in removing the 
2175A dip. 

As a consequence, the coagulation scenario is also sup- 
ported by the lack of a prominent carbon dip in spectra of 
reddened Syls (paper I). However, the depletion of small 
grains can also be achieved with other physical processes 
(eg. sublimation, Sect. 2. 5) and, therefore, the absence of 
the carbon dip alone cannot be used to rule out other 
scenarios. Yet, whatever is the mechanism responsible for 
depleting small grains, the absence of the carbon dip cer- 
tainly indicates that the properties of dust in AGNs are 
different from Galactic and, in particular, that the dust 
distribution is biased in favor of large grains. 

4. Conclusions 

We have discussed various interpretations of the observa- 
tional evidences for "anomalous" properties of dust in the 
circumnuclear region of AGNs and, in particular, the re- 
duced Eb-v/Nh, the reduced Ay/Nn, the absence of the 
silicate absorption feature in the mid-IR spectra of Sy2s 
and the absence of the carbon dip in the UV spectra of 
reddened Syls, which are reported in a companion paper 
(Maiolino et al. 2000, paper I). 

A dust distribution dominated by large grains as a con- 
sequence of coagulation is the interpretation which best 
fit the various observational constraints. The formation of 
large grains is naturally expected in the high density envi- 
ronment characterizing the circumnuclear region of AGNs. 
A dust grain distribution biased in favor of large grains 
makes the extinction curve flatter than Galactic and fea- 
tureless (especially for what concerns the silicate mid-IR 
feature and the UV carbon dip). Coagulation also yields 
a reduced Eb-v/Nh and Ay/Nn- 

It should be noted that nearly all of the papers, pub- 
lished so far, dealing with the effects of dust in AGNs have 
adopted a standard Galactic extinction curve (essentially 
with the only exception of Laor & Draine 1993, who first 
proposed the large grain scenario). The assumption that 
the properties of dust in the extreme conditions of the 
circumnuclear region of AGNs are similar to the Galac- 
tic diffuse interstellar medium is somewhat naive, given 
that even in the denser clouds of our Galaxy the extinc- 
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tion curve is already different with respect to the diffuse 
medium. Probably, many of the past studies considering 
the effects of dust obscuration and scattering in AGNs 
should be revised at the light of the findings obtained in 
this paper. 

With regard to the scattered radiation (especially in 
obscured AGNs) it is worth noting that a dust distribu- 
tion dominated by large grains would yield a nearly gray 
reflection, thus mimicking the scattering by free electrons. 
Large grains produced by dust coagulation would have a 
reduced scattering efficiency with respect to the standard 
Galactic dust mixture. However, while in some cases coag- 
ulation can result in a reduction of the scattering efficiency 
by more than an order of magnitude, for some of the dust 
grain distributions producing a flat extinction and scatter- 
ing curve the reduction in the scattering efficincy is much 
lower (e.g. for the model discussed in 3.2.2 the scattering 
efficiency is reduced by a factor of ~2). Therefore, it is 
possible that some of the gray-reflected spectra found in 
some obscured AGNs, and ascribed to electron scattering, 
might actually be produced by large grains scattering. 

We have also examined other possible interpretations 
for the "anomalous" dust properties discussed in paper I, 
and more specifically: 

1) a metallicity higher than solar (as often found in the 
central region of AGNs) might affect the estimate of the 
equivalent N# from the photoelectric cutoff in the hard 
X-rays; 

2) the observed broad lines might be scattered by a re- 
flecting medium through a path of much lower gaseous 
column than the X-rays; 

3) reduction of the scattering component of the extinction 
curve as a consequence scattering of the radiation into our 
beam; 

4) the BLR might have a size comparable to that of the 
absorbing medium (at variance with what assumed by the 
standard unified model) and, therefore, suffer less absorp- 
tion with respect to the nuclear X-ray source; 

5) the hard X-rays might be partly absorbed by the same 
clouds which produce the broad lines; 

6) the sublimation of dust at the inner face of the obscur- 
ing torus would produce a gaseous dust-free region; 

7) the absorbing medium might simply be characterized 
by a low dust-to-gas ratio (in terms of mass). 

Although some of these cases might apply to some of 
the objects examined in paper I, none of them can fit all 
of the observational constraints for most of the objects, at 
variance with the large grains scenario. 
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Fig. 2. Expected variation of extinction and emission line 
properties as a function of the dust-to-gas ratio relative to 
the Galactic value. The predictions of the models, which 
assume a standard dust size distribution, are shown by the 
thick solid lines. The upper two panels show the expected 
variation of the Eb-v/Nh and Ay/Nn ratios. The line 
and shaded area marked by "AGN observed" are the mean 
and ±1<7 dispersion value found in the sample of paper I. 
The thin dashed vertical lines indicate the range of gas-to- 
dust ratios that matches the observed Eb-v/Nh for ^68% 
of the objects (±1<t). The thick dotted lines indicate the 
effect of an increased dust-to-gas ratio due to a metallic- 
ity three times solar. In the bottom panel "H/3(torus)" is 
the line luminosity produced within the dusty torus while 
"H/3(max)" is what one expects if all the ionizing pho- 
tons are absorbed by the gas (i.e. the dust-free case). The 
light horizontal line constrains the maximum line emission 
produced by the torus (in terms of H/3(torus)/H/3(max)) 
based on the maximum covering factor of the Narrow Line 
Region (~8% ). The thin solid line in the bottom panel 
shows, for comparison, the result obtained with Cloudy. 
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Fig. 3. Same as Fig. 2, but in this case the dust-to-gas 
ratio is constant while the maximum size of the grain dis- 
tribution a max varies along the abscissa. The exponent 
of the grains size distribution q and the minimum size 
a m in are kept constant to the standard values of 3.5 and 
0.005^m, respectively. In this figure we also show an ad- 
ditional panel, the top one, which show the variation of 
Rv = Av/Eb-v- The thick dotted line in the bottom 
panel shows the effect of a higher dust-to-gas ratio due to 
a metallicity three times solar. 
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of the grain distribution a max is kept constant and equal 
to 1/zm while the exponent q is varied along the abscissa. 
The thick dashed lines show the case in which 15% of 
the gaseous absorbing column (the outer part) contains 
Galactic standard dust. 
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Fig. 6. Depth of the silicate feature relative to the fea- 
tureless mid-IR absorption as measured by ratio (T9.7 — 
T 7.7)/ T 7.7- The upper horizontal line gives the Galactic 
value, while the lower horizontal line gives the upper limit 
obtained for Sy2s (the shaded area gives the spread of the 
upper limits). The upper dashed curve gives the ratio ex- 
pected by the model presented in Sect. 3. 2.1 (q=3.5) and 
is related to the upper abscissa scale, while the lower solid 
curves give the ratio expected by the models presented in 
Sects. 3. 2. 2 and 3.2.3 (a ma x = 1/im and a max = 10/xm) and 
are related to the lower abscissa scale. 
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Fig. 5. Same as Fig. 4, but with a max = 10/mi. The thick 
dashed lines show the case in which 3% of the gaseous 
absorbing column (the outer part) contains Galactic stan- 
dard dust. 
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Fig. 7. Comparison of the Galactic UV extinction curve 
with the extinction curve given by a dusty medium which 
is depleted of grains smaller than 0.03/^m. Note the ab- 
sence of the carbon dip in the latter case. 



